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by Henry J. Hettel, Roman Krawec, George M. Prok, and Clyde C. Swett 

Lewis Research Center 

SUMMARY 

The admixture of helium to a hydrogen plasma atmosphere has increased the radio- 
frequency power absorption in the ion-cyclotron mode. Measurements of power absorp- 
tion, plasma radius, electron density, and magnetic flux density at the peak power point 
were obtained as a function of plasma composition. A relation was derived involving 
these experimental parameters and the ion density. It is shown that the observed power 
enhancement can be explained if it is assumed that the metastable states of the helium 
atoms react with hydrogen atoms and molecules to increase the plasma ion density. 


INTRODUCTION 

For efficiency in heating protons by means of an ion-cyclotron wave, the proton 
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density should be on the order of 10 per cubic meter (refs. 1 to 4). The power which 

can be coupled into the wave is dependent not only on the ion density, but also on the 

effective proton plasma radius. Thus, for a given experimental configuration, increased 

power can be coupled into the wave if a way can be found to increase the plasma radius 

and/or the ion density (provided the ion density does not already exceed the optimum). 

On theoretical grounds, the admixture of helium to the hydrogen support gas 

promised to provide such a means via collisions of the second kind between hydrogen 

species and metastable helium atoms. The 19. 8 electron volts (3.17 aJ) of the 
3 

helium 2 metastable state is sufficient to ionize hydrogen atoms and dissociate hydro- 
gen molecules. Thus, the helium metastables produced by electron impact in the primary 
electron beam region (refs. 5 and 6), can diffuse across the confining magnetic -field 
lines and produce ionization beyond the confines of the primary beam, thereby increasing 
the proton plasma radius. Simultaneously, an increase in proton density will result, 



the degree of increase depending on the magnitude of the proton production rate 
ascribable to helium metastables relative to that ascribable to electron impact. A series 
of measurements were made in the ion cyclotron resonance apparatus (ref. 5) using 
various mixtures of hydrogen and helium gas. Variations were observed not only in the 
power absorbed, but also in plasma radius, electron density, and resonance frequency 
as the gas composition was changed. A theoretical model is set forth whereby the role 
of the metastable helium atom in increasing the proton density may be deduced. 


APPARATUS 


A schematic diagram of the experimental apparatus (ref. 5) is shown in figure 1. 



Figure 1. - Configuration of ion-cyclotron resonance apparatus 
2B (ICRA 2B). 


The vacuum chamber is a 2. 7-meter-long by 10-centimeter-diameter tube consisting of 
an 84 -centimeter -long aluminum oxide center section with stainless -steel end sections. 
The axial magnetic field has a 2 to 1 mirror ratio, can operate continuously, and may be 
varied in strength to 1 tesla (10 000 G) in the flat -field portion between the mirrors. 

The magnetic -field strength at the cathode location is about 15 percent of the flat -field 
value . 

A commercial grade of hydrogen gas or a previously prepared and analyzed mixture 
of commercial grade hydrogen and commercial grade helium was fed into one end of the 
vacuum chamber and pumped out the other end. The flow rate, and thus the pressure, 
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Figure 3. - Schematic plasma power curve. 


was regulated by means of a variable leak valve upstream of the vacuum chamber. The 
pressure was measured by means of a McLeod gage mounted about 0. 9 meter downstream 
of the center of the test chamber, which is about 2. 7 meters from the hydrogen inlet. 
Under these conditions, there was a twofold pressure drop from one end of the test sec- 
tion to the other. For all experiments reported herein, the pressure of the test gas was 

O 

held constant at 2.00±0.05 microns (0.267±0.007 N/m ) as measured by the McLeod 
gage. 

A sketch of the hot cathode electron source used is presented in figure 2. Details 
of construction and operation have been previously published (ref. 6). 


THEORY 

Wave Power Input Related to Experimental Parameters 

A schematic plot of radio -frequency (RF) power input to the plasma with variation 
in magnetic -flux density is shown in figure 3, where P r is the transmitter output power 
at resonance and P Qr is the off -resonance or no -wave -propagation power, so that their 
difference P w is the power being fed into the ion cyclotron wave. The base power P^ 
is the power dissipated in the electrical components when there is no plasma. It is 
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measured by running the apparatus in the same fashion as for a normal test run, except 
that no support gas is admitted to the test chamber so that no plasma can exist. 

We shall be concerned primarily with the peak power point. According to Stix 
(ref. 7), the quantity P^P^ (defined by him as W) is given by 


= QUV 


where Q is the usual figure of merit for a coil (the ratio of its inductive reactance to 
its resistance) and U and V are defined as follows: 

u, - 8 * r2 ^ 2 | 

ttX 3 /m - - f 


Ij (kR)K 1 (kR ) 

V = ± L. 

Ij(kR JKjOeR) tS 


The notation used in these and subsequent relations are defined in appendix A. The 
approximate dispersion relation for ion-cyclotron waves is given by Stix (ref. 7, eqs. 10 
and 26) as 


0 = 1 + 


Z 2 nX 2 ; 
2. 03 x10 18 A 


Inserting the values of U and V given by equations (2) and (3) yields for equation (1) 


P w _ 8QZR 2 0 2 L fs Ij(kR) K^kR^ 


The gross behavior of the Bessel function ratio I^/K^ may be extracted by writing them 
as the product of a deviation term 8 and their small argument approximations 


ii(4) =£, Kj(o = - 
1 2 1 I 
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for £ < 1. Thus, 6(kR) is defined by the relation 


I^kR) 

Kj(kR) 


k 2 R 2 

5(kR) 

2 


which, on noting that k = 2ir/\, yields for equation (5) 



8QZ L fg 



5(kR) 

6(kR c ) 


R 4 J2 2 


( 7 ) 


( 8 ) 


In arriving at equation (2), Stix specified that the wavelength in the plasma was 
equal to the geometrical coil wavelength. A more exact specification might be that the 
plasma wavelength be equal, not to the geometrical coil wavelength, but rather to the 
apparent coil wavelength as seen by the plasma. For a perfectly wound infinite -length 
coil, the geometrical and apparent coil wavelengths are the same, but, for a finite 
length coil with end effects, the two may differ. Indeed, a Fourier decomposition of 
the RF vacuum magnetic field under our coil shows significant contributions from wave- 
length components of 89, 44. 5, and 29.6 centimeters (ref. 5). Thus, we have no a priori 
knowledge even of the gross magnitude of the plasma wavelength, let alone of what small 
changes in apparent coil wavelength such coil loading factors as proton density level and 
radial density profile may engender. 

Our ignorance of the apparent coil wavelength may be obviated by letting the plasma 
tell us what is the coil wavelength it sees. This is readily done by use of the dispersion 
relation. Thus, if equation (4) is written for the proton-cyclotron wave by setting 
Z = A = 1, then 



i 3 / 2 « 3 


9 3/2 1s 3/2 

(1 - ST) (2.03X10 10 ) 


which yields, upon substitution into equation (8), 
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Grouping the constants permits the writing of 



n 


3 /2rV 


5(kR) 

2 3/2 S(kR ) 

(i - sr) c 


(ii) 


from which it is seen that, in view of the low order dependence of £2 on n, indicated by 
equation (4), an increase in the value of n and/or R by helium admixture should 
increase the power coupled into the plasma wave. 

For a helium -hydrogen mixture whose composition is given by x, the mole fraction 
helium, the quantity (P /P K ) may be conveniently compared to the similar quantity for 

W D x 

a pure hydrogen plasma (x = 0) by defining a ratio F such that 


F = 


<y p b> x . < p w> v 

<V p b>0 < P w>0 


( 12 ) 


The base power P^ is measured under no-gas conditions, and it is independent of x. 
Thus, F is simply the wave power ratio between the mixture and pure hydrogen. 

If coupling is negligible between the proton-cyclotron wave and the He + cyclotron 
wave, then equations (1) to (11) may be written not only for pure hydrogen supporting 
atmospheres, but also for hydrogen-helium mixtures as well. Hence, in view of 
equation (11), equation (12) becomes 


F = 




3/2 

-n*/ \ 6<k o R o>/ 6 < k o R c>/ 


(13) 


It will prove convenient to define a quantity a by 

A . 6(k x R x )/5(k x R c ) 
6(k Q R 0 )/5(k 0 R c ) 


(14) 


For x = 0, A is unity. Asa first order approximation, take it to be unity for all other 
x as well. The validity of this approximation is demonstrated in appendix B, wherein 
it is shown to have no significant effect on the results of this investigation. Thus, 
equation (13) becomes 


F 


meta 




(15) 
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Experimental measurements are acquirable for the quantities R x /Rqj ^ x > and so 
that prediction of F depends only on the prediction of n^/nQ. The functional x- 
dependence of n x /nQ is derived on the basis of enhanced hydrogen ionization by helium 
metastables. Insertion of the result into equation (15), together with the experimental 
values of £2 and R x /Rqj will provide F me £ as a function of x. The correctness of 
the metastable interaction assumption may then be assessed by comparing F me ^ with 

^exp" 


Proton Density and Helium Metastables 


In evaluating the effect of the helium metastables on the proton production and loss 
rates, one must use a model consistent with the known characteristics of the plasma 

being studied. The plasma of this experiment has a total neutral density of about 

1 9 

5x10 per cubic meter. From electron density measurements the degree of ionization 
is in the vicinity of 0. 5 percent, and the degree of dissociation is estimated from visual 
spectroscopic observations to be between 5 and 25 percent. The mean free path of the 
molecular species is about 20 centimeters, whereas the discharge tube radius is 
5 centimeters. Thus the following assumptions are consistent with the plasma of this 
experiment: 

(1) Collisions between electrons and ionized species are insignificant because of 
their low concentrations with respect to those of hydrogen and helium. 

(2) Wall losses of ions and excited plasma neutrals predominate over losses through 
plasma-phase reactions; hence, loss rates of ions and metastable atoms will be directly 
proportional to their respective concentrations. Thus, the reactions occurring in the 
plasma are restricted to those shown in table I and have the rate constants defined 
therein. A discussion of why these reactions were chosen is presented in appendix C. 


TABLE I. - PLASMA REACTIONS 


Production reaction 

H 2 -£i»2H 

H 2 + He m — <~2H + He 0 
H—^-*-H + + e 
H + He m — >-H + + e + He 0 


Cross 

Loss reactions 

Rate 

section 


constant 

a i 

h 2 

K i 

a 2 

H wall V2 H 2 

K 2 

a 3 

Hem -^r He ° 

K 3 

a 4 



a 5 
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The duration of plasma operation being on the order of minutes, a steady state is 
assumed, permitting rate balance equations to be written for all species. For the gas 
phase reaction 


A + B— X + Y 


The production rate per unit volume, may be written 


dn(3Q _ dn(Y) 
dt dt 


= n(A) n(B)a AB v AB 


( 16 ) 


(17) 


where n(A) and n(B) are the number densities of the reactants A and B, ct ab is the 
cross section for the reaction in question, and v AB is the velocity of closure of A and 
B. As written, equation (17) ignores the fact that cr AB is a function of v AB> and that, 
for an assembly of particles, the average value of the product °ab v AB s ^ ou ^ be used. 
For electron impact reactions, a AB varies as log v AB - Anticipating that v AB will be 
found to vary only a little, so that cr AB will vary even less, we shall consider all a AB 's 
to be constants. Under these circumstances, v AB then becomes the average value of 
the electron velocity distribution. 

The total plasma production rate for any species is given by the product of the 
plasma volume and the production rate per unit volume. Thus, for a uniform plasma of 
radius R and length l 


S!M = n(A)n(B)<7 ABW R h 


(17a) 


where dN(x)/dt is the total plasma production rate. If the plasma losses are to be all 
to the wall, then all species formed in the plasma must pass through the plasma boundary. 
Thus, the total plasma loss rate is the product of the plasma surface area and the loss 
rate per unit plasma surface area. Thus, 

. MM = K(x)n(x)2ffRZ (17b) 

dt 

where K(x) is the loss rate constant of species x. 

For the reactions of table I, equating the proton production and loss rates yields 


dN(H + ) 

dt 


= n e n(H)cr 3 v+ n(H)n(He m )cr 4 w = 2K 1 n(H + )/R 


(18) 


8 



where v is the average closure velocity of electrons with plasma neutrals, and w is 
correspondingly the average closure velocity between plasma neutrals. Similarly, for 
the helium metastables, 


dN(He ) 
dt 


= n e n(He)ff 5 v = 2Kgn(He m )/R 


Elimination of n(He m ) between equations (18) and (19) yields 


{ a Q v or-wovvR 1 

-A. + JL — 2 n(He)> 

2K 1 4K 1 K 3 J 


(19) 


( 20 ) 


A relation exists between n(H) and n(He) which stems from the composition of the 
gas feeding the plasma. This relation may be expressed as follows: The total number 
density n(tot) of all constituents in the feed gas is given by 

n(tot)° = n(H 2 )° + n(He)° (22) 

where the superscript zeros denote feed gas to distinguish them from the corresponding 
quantities in the plasma which carry no superscript. In terms of the mole fraction 
helium x 

n(He)° = xn(tot) 0 (23a) 

n(H 2 )° = (1 - x)n(tot) 0 (23b) 

Furthermore, a quantity m, the degree of dissociation, may be defined as that fraction 
of the original molecular hydrogen in the mixture that has dissociated into atoms. Thus, 
since every hydrogen molecule yields two atoms on dissociation, 

n(H) = 2mn(H 2 )^ = 2m (1 - x)n(tot)® (24a) 


and 


n(H 2 ) = (1 - m)n(H 2 )° = (1 - m)(l - x)n(tot) 0 (24b) 

Having now obtained expressions for the plasma quantities n(H) and n(H 2 ), it re- 
mains only to express n(He). Since the loss of ground-state helium atoms due to ioniza- 
tion, metastable formation, etc. , is insignificant in the slightly ionized plasma of this 
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experiment, the helium density in the plasma is different from the helium density in the 
feed gas only because of the slight amount of neutral heating which occurs in the plasma. 
The temperature rise of the plasma neutrals is estimated to be about 50° C. Thus, 

n(He) = tn(He)° (25) 


where the factor t takes care of the reduction in plasma helium density resulting from 
thermal expansion of the plasma atmosphere. Equation (23a) then becomes 

n(He) = xtn(tot)® (26) 


Now that the relations between n(H), n(H 2 ), and n(He) resulting from feed gas 
composition are established, they may now be inserted into equation (20) yielding, after 
factoring the quantity agv/2Kj, 


n(H + ) 


Ogn(tot) 1 

K, 


Rvm(l - x)n g 


o\.ovn(tot)0wt 

1 + xR 

2cr 3 K 3 


(27) 


Of the quantities in equation (27), v, m, and n g vary significantly with x. The 
remaining quantities Oj, Kj, w, and t are, at best, second order in x, a fortunate 
happenstance because there is no experimental data on the variation of these with x. 
Thus, the following constants are defined as 


a = 


agn(tot)® 


K 


1 


(28a) 


r s 


a^agn(tot)°wtRQ 


2 ° 3 K 3 


(28b) 


with the aid of which equation (27) becomes 


R 

n x< H+ > ' “ v x m x (1 - + r jp- x ' 


(29) 


in which the x -dependence of the pertinent quantities has been indicated by subscripting. 
For a pure hydrogen plasma, for which x = 0, equation (29) becomes 



Thus, from equations (29) and (30), the ratio n x /n^ is 


n (H + ) v m n 


= _x _ex (1 _ x) _£L /1 + yx 


R_ 


n 0 (H + ) v 0 m 0 n e0 


R 


0 


R r 


( 31 ) 


An expression for m /m n is readily derived by considering the probable dissoci- 

X u 

ation reactions in the plasma, as discussed in appendix C. The rate balance equations 
for the atomic hydrogen production and loss reactions of table I yield the following 
relation: 


dN(H) 


dt 


= 2 a-jvngnfag) + 2cr2wn(H2)n(He m ) = 2K2n(H)/R 


(32) 


After rearrangement of equation (32) and use of equation (19) to eliminate n(He m ), 


n(H) = 2Rn e n(H 2 ) 


CT 1 V OVjWOVVR 

1 + J _5 n(He) 


2K 2 4K 2 K 3 


(33) 


When equations (24a), (24b), and (26) are used, equation (33) becomes 

v0 


™ °i v 

-2- = - 1 - n 


1 - m 2 Kr 


a 0 ov wtn(tot) 

1 + X R 

2aiK 3 


(34) 


The coefficient of x may be expressed in terms of y, defined by equation (28b) by 
writing 


a 2°5 W t n (t°t)^Ro i 

,a _2 fj\ 

cr 4 (J 5 W t n (t 0 t^Ro 


2 °-l K 3 ' 

"4/ 

2cr 3 K 3 



•y 


(35) 


Thus, by defining the quantities 


^2 ^3 

a l °4 


cr. 


2K 0 


(36a) 

(36b) 
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and using x-subscripts to designate the x-dependent parameters, equation (34) may be 
written as 


m 


R_ 




Writing equation (37) for the pure hydrogen case yields 


m. 


= |v 0 ne 0 R 0 

1 - m 0 


( 37 ) 


(38) 


whereupon the ratio of equations (37) and (38) is 



(39) 


Thus, on substitution of equation (31) into equation (15), F met is given by 



R 

(1 - x) (l + — yx) 
Rn 


13/2 


(40) 


EXPERIMENTAL RESULTS 
Wave Power Measurements 

A composite plot of total input power against magnetic flux density in the vicinity 
of the peak power point is presented in figure 4 for hydrogen and the six mixtures 
studied. Gas composition is indicated by x, the mole fraction of helium in the mixture. 

No variation in plasma power input with magnetic flux density was observed for 
the pure helium plasma x = 1 over the 0. 3 to 0. 5 tesla range studied. This justifies 
the assignment of the plasma power input at 0. 32 tesla 1 in the mixture under study as 

1 No significance should be attached to the specific P or value at 0. 32 tesla, since 

the P is constant over the range from 0. 30 to 0. 36 tesla. The 0.32 -tesla value 
or 

happened to be convenient to establish in our device. 


12 




Magnetic flux density, T 
Figure 4. - Experimental plasma power curves. 

being also the off -resonant power input P Qr at the peak power field for that mixture. 
Thus, the wave power input at any flux density B is given by the difference between the 
total-power input at that flux density P(B) and the total power input at 0. 32 tesla 
P(0. 32); that is, 

P w (B) = P(B) - P(0. 32) (41) 

Hence, for B = B (peak) of each mixture 

p w = P r " P(0 * 32) (42) 

For pure hydrogen, P w was 950 watts. 

The quantity F , defined by equation (12), is plotted against gas composition in 

0Xp 

figure 5. The enhancement of wave power input by helium admixture is seen to peak 
near x= 0.6. 
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Measurement of Coil Driving Frequency to Ion Cyclotron Frequency Ratio 


The quantity £2 is defined by Stix (ref. 7) as the ratio of the angular frequency of 
the propagated wave to the angular ion-cyclotron frequency, the latter being given by 
ZeB/nip. For the experimental case where the driving frequency co is fixed, £2 becomes 
the ratio B c ^/B, where B c - = unip/ Ze and B is the magnetic flux density at which 
maximum power transfer to the wave occurs. The values of £2 so obtained are plotted 
against x in figure 6. 

Another, and independent, measurement of £2 was also made, however. A magnetic 
probe sensitive to the axial component of the plasma wave was positioned downstream of 
the RF coil. The construction, operation, and positioning of this probe is described 
in reference 5. The amplitude of the signal from this probe is a measure of the strength 
of the wave passing through the plasma, and hence the peak probe amplitude is also an 
indicator of £2. The values of £2 obtained from the probe measurements are also 
plotted in figure 6. 

Except for the x = 0. 22 and x = 0. 328 points, the agreement between the two 
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measurements of £2 is good. Since we believe that the probe is a more reliable indicator 
of O than the power measurements, we shall in all cases use the probe values in the 
calculations. 


Plasma Radius Measurements 

The radial distribution of the radiation emitted by the plasma was obtained by the 
standard technique of performing an Abel inversion on the measured chord intensity 
distribution (ref. 8). The device for obtaining the measured chord intensity distribution 
(fig. 7) is basically a two-slit collimator with the slits parallel to the plasma axis. The 
cross-sectional dimensions of the plasma slab viewed are 0. 75 by 10 millimeters. Since 
no dispersing element was interjected between the photocathode and the collimator of the 
scanner and since the viewing port window was fused quartz, the photomultiplier saw all 
radiation within the envelope of its spectral sensitivity curve. 

A plot of the radial intensity distribution of the plasma is shown in figure 8 for the 
peak power point of the x = 0. 50 mixture. This distribution is typical of all those 
observed in this study. The significance of the intensity distribution is that, in a tenuous 
and slightly ionized plasma such as ours, the radiant intensity emitted from a point is 
directly proportional to the electron density at that point. 

Two problems exist in interpreting the radial distribution curves: (1) The quantity 
of interest is the effective proton plasma radius, whereas the quantity measured is the 
radial distribution of electrons possessing sufficient energy to excite hydrogen molecules 
and helium atoms, and (2) the radial distribution is a smooth function for which the 
designation of a single value characterizing an effective radius is not free of ambiguity. 
With respect to the first problem, the assumption will be made that the effective proton 
plasma radius is directly proportional to the effective electron plasma radius since, 
after all, the ions are generated by ionization of neutrals by the same energetic plasma 
electrons that are responsible for the excitation of the neutrals. With respect to the 
second problem, the intensity-weighted average radius of the distribution R, given by 



is used as the characteristic radius of the plasma. The quantity I(r) is, of course, the 
radial intensity distribution of the plasma as derived from the Abel inversion of the 
chord intensity distribution. 

The quantity to be inserted into equation (15) is the proton plasma radius ratio 
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Radiant intensity arbitrary units 



Figure 7. * Schematic diagram of radiation scanner. 



0 1 2 3 4 5 

Plasma radius, cm 

Figure 8. -Typical experimental radial intensity distribu - 





R x /R 0 . If R is proportional to the effective electron plasma radius, and if this, in 
turn, is proportional to the proton plasma radius, then the error in R x /Rq will be 
less than the error in the absolute value of the effective plasma radius. To illustrate 
this, two other methods of obtaining a characteristic radius were investigated: (1) The 

p 

calculation of R in the standard manner analogous to equation (43), and (2) the calcu- 
lation of an equivalent flat profile (dashed line in fig. 8) such that area A equals area C 
equals one -half area B. Although the magnitude of the three characteristic radii 
differed greatly for a given radial distribution, curves of R x /Rq calculated from each 
technique were coincident with one another to within a few percent. 

A plot of R x /R Q at the peak power point, derived from the R values is presented 
in figure 9. 



Figure 9. - Experimental radius ratios. 



Electron Density Measurements 


An 8-millimeter microwave interferometer was used to measure space average 
values of the electron density. Complete details of the system and its operation are 
given in reference 6. The plot of the electron density ratio at the peak power point is 
given as a function of gas composition in figure 10. 


Electron Velocity Measurements 

The total potential drop required to sustain a 15 -ampere beam current was another 
of the experimental parameters measured. It was shown in reference 6 that the electron 
energy in the primary beam was directly proportional to the beam power input. Thus, 
since the current was held constant for all runs, the electron energy is proportional to 
the measured potential drop. It follows that 
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( 44 ) 


where E is the potential drop across the plasma whose composition is characterized 
by its mole fraction x. The electron velocity ratio at the peak power point is given as a 
function of gas composition in figure 11. 



Figure 11. - Experimental electron velocity ratios. 


DISCUSSION 


Equations (39) and (40) can, in principle, be used to calculate F 


met 

the gas mixtures studied. Experimental values are available for £2 X , £2 q, 
ratios R x /R q , n ex /n e0 , 


for each of 
and the 

and v x /vq. Thus, if mQ, /3, and y were known, m x could 
be calculated from equation (39) and then F me ^. obtained from equation (40). 

Since nip, the degree of dissociation of the pure hydrogen plasma, is not known, 
the range of its probable value must be determined. By definition mQ is a number 
between zero and unity. The best estimate for its value in the plasma of this experiment 
is 0. 15, so its probable range is between 0 and 0. 5. 

A discussion of the probable range of the cross-section parameter j3 is presented 
in appendix C, where it is shown that the most probable value of j3 lies within the range 
0. 01 s p < 100. 

Because so little is known about the quantities comprising y, its probable magnitude 
cannot be estimated. No recourse exists, therefore, but to determine y by equating 

The experimental point x = 0. 805 was arbitrarily chosen. 

x’ S’ R x /R 0’ n ex /e 0’ and v x /v 0 
used were the raw experimental values presented in figures 6 and 9 to 11, respectively. 

To eliminate division-by-zero problems in the computer program written to calculate 


F’exp at one point. 

In the calculations to follow, the values of Q 


F meta to 
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(a) (3 = 3. (b) (3 = 10. 

Figure 12. - Comparison of calculated and experimental wave power ratios for various values of m Q and (5. 

F met’ the condition m Q = 0 was approximated by setting m Q = lxlO -6 . Figure 12(a) 

portrays the results of setting /3 = 3, deriving y by equating F , to F at the 

m g t g xp 

X = 0. 805 point, and letting m Q take the values 0, 0. 01, . . . , 0. 5 to yield the family 
of curves shown. For this value of / 3 , F met does not approach F except at the 
match point. This is true for all smaller values of /3 as well. 

Figures 12(b) to (d) demonstrate the effect of increasing /3 to 10, 30, and 100, 
respectively. As /3 is increased, concordance of F met and F exp occurs at continu- 
ously lower values of niQ, F gx p being close to the estimated value for the plasma when 
/3 is in excess of 30. 

As different values of the constants m Q and /3 are inserted into equations (39) and 

(40), different values of the constant y are also required to equate F met to F 

at x = 0. 805. These values of y are reported in figure 13. There is marked on each 

curve the y, m n pair which yields best concordance between F m . and F for the 
u met exp 

value of /3 chosen. Thus, the line joining these points is the locus of the (3, y , m Q 
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(c) (3 = 30. (d) 0 = 100. 


Figure 12. -Concluded. 


solution set for the model used to derive F me ^. . This locus does not depend on which 
experimental point is chosen for the gamma -determining fit, because the best- 
concordance F me j. curves agree well with F gX p at all points below x = 0. 805. The 
agreement at x = 0. 959 is poor, probably because the model is too simple to adequately 
treat the case of such small hydrogen concentrations. 

In this proposed model, no other effect has been ascribed to helium save the 
interaction of its metastables with molecular and atomic hydrogen. Hence, the ability 
of this single assumption to explain the gross features of the experimental power en- 
hancement lends support to the belief that this mechanism plays a major role in the 
phenomenon observed. 

The addition of helium to hydrogen will not result in increased RF power coupling 
to the plasma in all cases, however. The metastables enhance power addition by provid- 
ing both increased proton density and increased plasma radius. Thus, it follows that 
no enhancement of power will occur on addition of helium to a hydrogen plasma whose 
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Figure 13. - Relation of y parameter to m c and |3. 

density is already optimum and whose radius is already the maximum permitted by the 
discharge container. 

That the helium does not act simply as a neutral diluent is readily demonstrable. 
Were it to do so, the F -behavior of the mixtures would be described by equations (39) 
and (40) with y = 0 and R x /Rq = 1 . Setting y = 0 erases only the effect of the proton 
density enhancement by the helium metastables. This results in the curves labeled 
F ffr = 0) figure 15. When, in addition, the increased radius effect produced by 
the metastables is also erased by setting R x /Rq = the curves labeled F me j. 

(y = 0, R X /R Q = 1) result. This second set of curves clearly disagree with the experi- 
mental measurement. 

Admittedly, it has not been unequivocally proven that the helium metastables are 
responsible for the enhanced power absorption in hydrogen-helium mixtures. Admittedly 
too, the evidence would have been stronger had it been possible independently to arrive 
at the values of (3 and y. Nevertheless, it appears that an explanation for the power 
enhancement invoking an increased proton concentration through helium metastable 
interaction is a reasonable explanation of the experimental facts. 
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SLOT DISCHARGE 


As discussed in reference 5, the RF electrostatic shield used in this apparatus was 
rolled from copper sheet leaving about a 1. 8-centimeter slot the full length of the shield. 
It appeared that a slot of this width permitted some leakage of the E„ component of the 
electric field of the RF coil to occur, since an enhancement of the radiant plasma inten- 
sity is observed in the immediate vicinity of the slot. This so-called slot discharge may 
be the origin of the small, but finite, difference between the off-resonance plasma 
power, P Qr in figure 3, and the base plasma power which is coupled into the plasma 
The plasma radius measurements were made on the plasma hemicylinder opposite the 
slot discharge. 

During the preparation of this report, an improved shield design was installed which 
eliminated the slot discharge and reduced the quantity P Qr - P^ to zero. Enhancement 
of power input by the use of hydrogen-helium mixtures occurred nevertheless, indicating 
that the slot discharge played no dominant role in the enhancement process. 


CONCLUDING REMARKS 

The ability to couple radio frequency power into a hydrogen plasma is dependent 
upon many parameters, two of which are proton density and plasma radius. In studying 
the effect of support gas composition on these two parameters, it was noted that 

1. The admixture of helium to the hydrogen plasma atmosphere permitted more 
power to be coupled into the proton-cyclotron wave. The helium mole fraction of the 
optimum gas composition was 0. 65, for which the wave power input was 3.2 times as 
great as that for hydrogen alone. 

2. An almost linear increase in plasma radius with mole fraction helium x was 
noted over the range 0 < x < 0. 7, with the radius being observed to maximize near 
x = 0. 8. 

3. The assumption is made that. helium atoms in metastable states enhance the 
proton density in the plasma by collisions of the second kind with hydrogen atoms and 
molecules. On this basis, the dependence of the power enhancement on helium content 
was calculated, and it was found that this mechanism could account for the enhancement 
experimentally observed. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 11, 1967, 

129-01-05-09-22. 
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APPENDIX A 


SYMBOLS 

(The rationalized MKS system of units is used throughout unless otherwise noted.) 

A mass number of ion being resonated, in amu 
B magnetic flux density 

B z component of magnetic flux density along the z-axis 
E dc potential drop driving the primary beam 
e charge of electron 

e* excitation by electrons 

F eX p experimental wave power ratio using experimental power values 

F me t calculated wave power ratio based on metastable interaction mechanism defined 
by equation (15) 

Ip Bessel functions in the notation of Watson (ref. 9) 

I(r) radial intensity distribution of the plasma 

I(Y) ionization potential of species Y 

K loss rate constant (see reactions of table I) 
k wave number, equal to 2it/X 

Lj g ratio of the inductance of the fundamental component of the coil current to the 
total inductance of the coil system 

l coil length 

M mode number of the plasma oscillation 

m degree of dissociation, defined by equation (24a) 
m mass of proton 

r 

N total number of indicated species in the plasma volume 
n proton density, when species is unidentified 

n(4) number density of species 

P transmitter output power 

P^ no-plasma base power 
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pQ r off-resonant plasma power 

P r transmitter output power at resonance 
P w plasma wave power 

Q the ratio of inductive reactance of coil to its resistance 
R plasma radius 

R defined by equation (43) 

R coil radius 

t defined by equation (25) 

U defined by equation (2) 

V defined by equation (3) 

3 

V m energy of the helium 2 S metastable state 
v closure velocity of electrons with plasma species 

w closure velocity of helium metastables with plasma species 
x mole fraction helium 

Z charge number of ion in ion cyclotron wave 
a defined by equation (28a) 

j3 defined by equation (36a) 

A defined by equation (14) 

£ defined by equation (36b) 

6(kR) defined by equation (7) 
y defined by equation (28b) 

X wavelength of propagated wave 

ratio of coil driving frequency to ion cyclotron frequency 
co transmitter frequency, 41x10° rad per sec 

a reaction cross section 

£ dummy argument 

Subscripts: 
ci ion cyclotron 

e electron 
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j dummy, designating any member of the class under consideration 
m metastable state 

x the gas mixture whole mole fraction helium is x 

0 the gas mixture whose mole fraction helium is 0, that is, pure hydrogen 

Superscripts: 

* excitation by electrons 
m atom in m eta stable electronic state 
o atom in ground electronic state 
0 feed gas 
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APPENDIX B 


EVALUATION OF A 

If the definition of 6(kR), given by 


6(kR) = 


2 IjtkR) 
A^KjfkR) 


is substituted into the definition of A, given by 

A _ 6(k A )/5(k x R c ) 

6(k 0 R 0 )/6(k 0 R c ) 


( 7 ) 


(14) 


then 


M 2 ■iWiW K iW K M> 
V 'iWiW K iW K iW 


(Bl) 


To evaluate the Bessel functions, absolute numerical values must be obtained for the 
arguments. From the dispersion relation (eq. (4)), 



n 


, 1/2 


2. 03x 10 18 (1 - ft 2 ) 


(B2) 


so that, from the definition k = 2ir/\, 


k x = 2 ttO x| 


-il/2 


n 


2.03x10 


18 


(' - n xl 


(B3) 


17 

Absolute values of n x are not known, but nQ is known to be near 1x10 per cubic 
meter. Hence, if the relation between n x and nQ is known, we may write 


k x = 4.41X10" 9 n Q 


1/2 


K /n oK 


i - n! 


nl/2 


(B4) 
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TABLE n. - VALUES OF A 

Mole 
fraction 
helium, 
x 

0 

. 220 
. 328 
.495 
.580 
. 805 
.959 

where n x /riQ will be taken to be given by equation (31). 

Absolute values for the proton plasma radius must also be specified, which are most 
conveniently written in the form 

R x = R 0 (R X /R 0 ) (B5) 

R x /R 0 being given in figure 9. No alternative exists but to use R for the proton plasma 
radius, so that Rq = 1.03 centimeters. The mean coil radius R c is 7. 3 centimeters. 

All quantities being now specified, A = A(x) may be calculated. The results are 
given in table n, which show that A(x) does not deviate appreciably from unity. However, 
were the x -dependence of A taken into account, a closer fit of F me j. to F 0X p would 
result. Thus, the approximation A = 1 does not introduce any degree of conformity 
between theory and experiment. 


Wave 

number, 

V 

A 

4.386 

1.000 

4. 344 

1.003 

4.275 

.997 

3.639 

1.022 

3. 521 

1.035 

2.843 

1.051 

1.223 

1. 126 


27 




APPENDIX C 


PROBABLE HYDROGEN PLASMA REACTIONS 
Electron Impact Processes 

If possessed of sufficient energy, electrons colliding with hydrogen molecules may 
generate a variety of reactions. Consider the partial energy level diagram of the hydro- 
gen molecule (ref. 10) presented in figure 14. The spectroscopic notation used is defined 
in reference 11. Electron impact excitation may produce any of the excited states shown. 
The net effect of exciting these states are tabulated in table HI. 

The first excited state (in order of increasing energy and decreasing cross-section) 

3 + 

is the Isa 2pa 2 U state, which is dissociative. Thus, any molecule excited to this 
state will immediately dissociate into' two ground state atoms. Furthermore, the ex- 
citation of the molecule to any of the bonding triplet states will ultimately result in dis- 

3 3 

sociation also, because all and n states produced can radiate to the dis- 

S & 

q 3 3 

sociative 2 , and all bonding 2 U and n u states can radiatively cascade through 

the lsa2sa 2 state to the dissociative 2 . Thus, all triplet state excitation of 
g u. 



Figure 14. - Potential energy curves for some typical states of hydrogen molecule. 
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TABLE HI. - HYDROGEN PLASMA REACTIONS 


Reaction 

number 

Reaction 

Cross 

section 

le 

_* o Spontaneous 

H 2\ * ^( 1S a2sa 3 ^) dlssociation > 2H 

\ . Xi/ 

>Hg(bound triplets) ^ 

ct i 

lm 

He m n , Spontaneous 

H 2 -H*(lsa2sa 3 2 + u ) dissocial -- 2« 

ff 2 

2e 

H 2-^ H 2 + e 

°6 

3e 

H 2 -£^-H + H + + e 

°7 

4e 

H »-H + + e 

°3 

4m 

H Ji^H + + e 

°4 


molecular hydrogen will ultimately produce a net dissociative reaction. In table III, 
the short hand description of this process is designated as reaction le. 

Electron excitation of ground state molecules to excited singlet states can occur 
too, but these can all return radiatively to the ground state either directly or through 
a cascade path. Thus, excitation of singlet states results only in radiative energy 
losses, but no dissociation. There being no useful result, the singlet state excitation 
process is not depicted in table HI. 

The preceding discussions are predicated on the premise that electrons are colliding 
with hydrogen molecules not only in their ground electronic state, but in their ground 
vibrational state as well. Under conditions of adiabatic excitation, that is, where the 
excitation occurs so quickly that the molecular nuclei are unable to move much in the 
elapsed time, the internuclear separation of the excited state produced will be closely 
the same as that of the target molecule at the instant that excitation occurred. If 
attention is now directed to any of the excited state Morse curves in figure 14 and 
especially to the ground state of the hydrogen molecular ion (Hg), it may be observed 
that excitation at sufficiently small internuclear separation could yield excited states 
with sufficient vibrational energy to be dissociative. However, as a result of a program 
undertaken to calculate the electron impact cross sections of the excited states of 
molecular hydrogen (ref. 11), it was found that the fraction of ground state molecules 
having internuclear separations small enough to produce dissociative excited states 
was negligible. 

Applying these considerations to ionization of the hydrogen molecule leads directly 
to the conclusion that the products of the ionizing electron impacts will be almost wholly 
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Hgj with very little H and H + , the products of a dissociative ionization. Thus reaction 
2e of table m should have a high cross section and will thus be important, whereas 
reaction 3e should have a low cross section and thus be unimportant. Indeed, the 
experiments of Tate and Smith and of Newhall, as reported by Goodyear and Von Engel 
(ref. 12), bear out these conclusions. 

Having surveyed the possible reactions of molecular excitation and failing to find 
a proton source among them, it follows that the proton source must be ionization of the 
atomic hydrogen produced by reaction le. Thus, reaction 4e of table m is assigned as 
the proton source in a hydrogen plasma, not only by elimination, but also because it is 
a high cross section reaction, to which reference 11 testifies. 


Helium Metastable Interactions 

The helium 2^Sj and 2 *Sq metastable states have energies of 19. 8 and 20. 6 eV 
(3. 17 and 3. 30 aJ), respectively. Thus, more than enough energy is available to permit 
these metastables to dissociate hydrogen molecules or ionize hydrogen atoms on col- 
lision with them. Considering the second possibility first, examine the reaction 

He m (2 3 S) + H — He 0 + H + + e + 6. 3 eV (Cl) 

and enter it as reaction 4m in table HI. The cross section for this reaction is not 
known. However, Penning reactions such as this have been known for many years, and 
it will be instructive to inquire into the magnitude of the cross sections of those reactions 
which have been studied. For the reaction 

He m (2 3 S) + x — He 0 + Y + + e + energy excess (C2) 

the energy excess is given by 


Energy excess = V m - I(Y) (C3) 

where V is the potential energy of the metastable state and I(Y) is the ionization 

rn |-yj 

potential of the species Y. Figure 15 is a plot of the cross section for the He - Y 
reaction against the energy excess of the reaction, the data being taken from refer- 
ences 13 and 14. One observes that the cross section for the reaction appears to 
increase linearly with energy excess. Since there is no reason to expect that atomic 
hydrogen should behave any differently from argon, krypton, or xenon as a Penning 
gas, then the cross section of reaction 4m might be predicted to be about 10x10 ^ 
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Energy, aj 

Figure 15. - Cross-sections of various gases for Penning 
ionization by helium triplet metastables. 


square centimeter since its energy excess is 6. 3 eV (1. 01 aJ). This is a sizable cross 
section and, if true, would support the assumption that Penning effect ionization of 
atomic hydrogen is a significant source of protons in our plasma. 

The Penning ionization is, of course, simply a specific case of a general class of 
reactions, namely, collisions of the second kind, where the potential energy of the 
incident particle is transformed in significant measure to potential energy of the target 
particle, with the excess providing kinetic energy of separation of the reaction products. 
Thus, when a metastable collides with a molecule, it is possible that electronic ex- 
citation of the molecule may occur, not only into the ionization continuum, but also into 
an excited repulsive state from which dissociation will result. Groth and Oldenberg 
(ref. 15) propose this mechanism to explain the enhancement of nitrogen dissociation by 
excited krypton states. 

Furthermore, Beutler and Eisenschimmel (ref. 16) showed that neon metastables 
produced dissociation of molecular hydrogen. Thus, for the interaction of helium meta- 
stables with molecular hydrogen the analog to reaction le is added to table ni as re- 
action lm. From figure 12(d), the excitation potential of the repulsive ^2* state at 
0. 074 nanometer is about 11 eV. 
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Unfortunately, no quantitative data is available for dissociative reactions between 
metastables and molecules. 


Estimation of (3 

Fortunately, it is not necessary to have numerical values for the cross sections of 
table III. What is necessary is only to estimate the value of /3, which is defined by 
equation (36a) to be 



a l 


(C4) 


Thus, /3 is the product of two ratios. To evaluate it, consider first the ratio o^/ cr^ , 
which involves the cross section for the reactions lm and 4m of table in. The Penning 
ionization reaction (shown as eq. (Cl)) involves basically the formation of a collision 
complex between two particles which subsequently decomposes into three particles, 
namely, a helium atom, a proton, and an electron. But reaction lm involves also the 
formation of a collision complex between two particles which subsequently decomposes 
into three particles, namely, a helium atom and two hydrogen atoms. Thus, both re- 
actions involve the transformation of the potential energy of the metastable state into 
potential energy required to break bonds in the target species plus kinetic energy to 
separate the product species. In this light, reactions lm and 4m are similar, whereupon 
the assumption that a 2 is within an order of magnitude of seems valid. Thus, 

0. 1 < a 2 /cr 4 < 10 (C5) 

Regarding the ratio 03 /ap data for both these cross sections as a function of mono- 
energetic electron energy exist in the literature (refs. 17 and 18). After these data have 
been transformed to yield the cross sections as a function of Maxwell averaged electron 
energy, the electron energy dependence of the ratio crg/cr^ may be found. The ratio is 
0. 01 at 1 eV (0. 16 aJ) and increases monotonically to 8 at 50 eV (8 aJ). Because this is 
the energy range that has been measured in the plasma by probe and spectrometric tech- 
niques (ref. 6 ), it is assumed that 


0. 01 < ag/CTj < 8 (C6) 

Thus, we shall expect to find /3 within the range 

0.001 <j3< 80 (C7) 
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